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Abstract 

Mathematical methods previously used (Phillies, J. Chem. Phys., 122 224905 (2005)) to interpret 
quasielastic light scattering spectroscopy (QELSS) spectra are here applied to relate diffusing wave 
spectroscopy (DWS) spectra to the moments X 2n of particle displacements in the solution under 
study. DWS spectra of optical probes are like QELSS spectra in that in general they are not 
determined solely by the second moment X 2 . In each case, the relationship between the spectrum 
and the particle motions arises from the field correlation function gg (t) for a single quasi-elastic 
scattering event. In most physically interesting cases, <?« (t) receives except at the shortest times 
large contributions from higher moments X(t) 2n , n > 1. As has long been known, the idealized 
form g^\t) = ex.p(—2q 2 X(t) 2 ), sometimes invoked to interpret DWS and QELSS spectra, only 
refers to (adequately) monodisperse, noninteracting, probes in purely Newtonian liquids and is 
erroneous for polydisperse particles, interacting particles, or particles in viscoelastic complex fluids. 
Furthermore, in DWS experiments fluctuations (for multiple scattering paths of fixed length) in 
the number of scattering events and the total-square scattering vector significantly modify the 
spectrum. 
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I. INTRODUCTION 



Three decades ago, Hallett and students^ used quasielastic light scattering spectroscopy 
(QELSS) to study dilute, intensely scattering probe spheres diffusing through polymer so- 
lutions. The spheres were used as optical probes of the polymer solution's rheology. By 
employing the Stokes-Einstein equation 

d = M- m 

P 

an apparent (micro) viscosity rj^ was inferred from the diffusion coefficient D p . Here &b, T, 
and R are Boltzmann's constant, the absolute temperature, and the probe radius, respec- 
tively Hallet, et al.'s rj^ were nearly equal to the viscosities 77 determined with classical 
instruments. These experiments were a logical continuation of earlier work of Laurent, et 
al.— , who used ultracentrifugation to study the sedimentation of mesoscopic colloids through 
various complex fluids. 

QELSS, fluorescence recovery after photobleaching, and related techniques have since 
been used to study the diffusion of optical probes through a wide variety of complex fluids, 
such as polymer^, surfactant, colloid^, and protein^ solutions. These studies are united 
under the cognomen Optical Probe Diffusion (OPD). Effective ways to control the binding 
of at least some polymers to some probes^ are known, permitting calibrated measurements 
of the microviscosity and the measurement of the thickness of the bound polymep^. It 
was early recognized that the microviscosity and the viscosity 77 obtained from capillary 
viscometers are often not the same. In a few cases, exotic reentrant behavior is observed^, 
in which rj^ 7^ 77, but only over a narrow band of concentrations. Systematic studies of 77^, 
as a function of the concentration and molecular weight of matrix polymers in solutio n 12 ' 1 ? , 
investigated whether optical probe diffusion is only an expensive replacement (albeit helpful 
with microscale samples) for an inexpensive capillary viscometer, or whether it gives inter- 
esting physical results. Fortunately, the latter case is correct: Optical Probe Diffusion gives 
novel information on polymer dynamics. 

Artifacts such as probe aggregation or polymer adsorption by the probe cause the probe to 
move too slowly, leading to 77^ > 77, thus permitting a ready separation of several important 
classes of artifact from physically interesting results. In the absence of artifacts, it is generally 
the case that ry M = 77 or that ry M < 77 (sometimes 77 M 77). It is incorrect to propose that 
77 > r)n could arise from the non-zero shear rate of capillary viscometers. For 77 > i]^ to 



arise from viscometer shear, the polymer solution would need to be shear thickening. Shear 
thickening is a rare phenomenon not encountered with most polymers that have been studied 
with probe diffusion. 

QELSS measures the field correlation function gi 1 \q, t) of the light single-scattered by the 
probes. QELSS spectra have recently sometimes been interpreted under the cognomen mi- 
crorheology, in which it is claimed that QELSS spectra are related to particle displacements 
via 

gi 1 \q,t)=eM-q 2 XWm, (2) 
q being the scattering vector determined by the laser wavelength and apparatus geometry. 
In microrheology, from X(t) 2 , a time- dependent diffusion coefficient is inferred from the 
supposed form X(t) 2 = 2tD(t). Finally, a time- or frequency-dependent generalization of 
the Stokes-Einstein equation and time- frequency transforms are used to infer from D(t) 
a frequency- dependent viscosity, or storage and loss moduli, corresponding to the original 
DWS spectrum. This process for determining G'{uS) and G"{uj) from X(t) 2 may be compared 
with determinations using classical mechanical means. 

Equations superficially similar to eq El appear in Berne and Pecorai*. However, as is 
made clear in Berne and Pecora, equation El only refers to the diffusion of identical, true 
Brownian particles whose motion is governed by the Langevin model. In this model, there 
are no memory effects such as those arising from solvent viscoelasticity. For Brownian 
particles that satisfy the Langevin model, X(t) 2 = 2Dt, where D is the time-independent 
diffusion constant. In the QELSS literature, eq El is only applied to certain simple cases, 
e.g., monodisperse polystyrene spheres in water. 

The motions of diffusing particles in viscoelastic fluids do not follow the Langevin equa- 
tion, because the random forces have non-zero correlation times. Correspondingly, it is 
incorrect to apply eqElto interpret QELSS spectra of probes in viscoelastic fluids. Equa- 
tion El also incorrectly predicts \og(gi 1 \q,t)) ~ q 2 as a uniform result. Modifications to 
the Langevin equation for physically realistic, i.e., viscoelastic at adequately high frequency, 
motions are encompassed by the Mori formalism^, which moves from the fundamental, mi- 
croscopic Liouville equation for particle motion to a reduced-variable description in which 
the random force and corresponding drag coefficient are described by memory functions. 

In order to clarify the interpretation of QELSS spectra, we recently examined^ the 
relationship between gi (q, t) and the mean-n th moments X(t) n of the particle displacement 
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along a single coordinate axis. The higher-order moments incorporate fluctuations in the 
particle displacements around their mean-square values. We^ showed: The odd moments of 
X vanish. In general, all even moments X(t) n of the displacement contribute to gi (q, t). A 
simple diagnostic identifies the special case in which eqElis correct: If the QELSS spectrum 
agrees with eqEl then from Doob's First Theorem^ the QELSS spectrum is necessarily a 
pure exponential with a time-independent diffusion constant. Doob's theorem is a purely 
mathematical result on zero-correlation-time random walks, and uses no physical arguments 
or rationales. Conversely, if the light scattering spectrum is not a single exponential, then the 
contrapositive form of Doob's theorem shows that egH cannot describe the QELSS spectrum; 
gs(q,t) contains non-trivial contributions from all higher-order even moments X(t) 2n . 

Equation El describes the motion of a solution of Brownian particles in an ideal simple 
fluid. If the particles are not monodisperse, or if the fluid is viscoelastic, the spectrum ceases 
to be a simple exponential, and gi (t) does not follow eq. It is not possible to mask the 
contributions of the higher moments of X by going to small q. While eq|2] cannot be used to 
interpret light scattering spectra of probe particles in most complex fluids, nothing is wrong 
with the underlying QELSS spectra that eqElhas been used to analyze. Reanalysis of those 
spectra using conventional means such as cumulants or mode decomposition may well lead 
to useful information. 

Recently, an alternative technique that uses light scattering to study the diffusion of 
optical probes has been propose d 18119 . Diffusing Wave Spectroscopy (DWS) studies light 
that has repeatedly been scattered by an intensely-scattering, colorless, opaque complex 
fluid. In DWS experiments considered here, scattering is due to diffusing mesoscopic particles 
that lead to multiple scattering of the incident light. The individual scattering events are 
describable as quasi-elastic single scattering; scattering events are far enough apart from 
each other that the motions of particles in a multiple scattering series are uncorrelated. 

The underlying experimental apparatuses in QELSS and DWS are the same, with a 
laser, scattering sample, and photomultiplier tube or equivalent photodetector. Differences 
appear in the scattering cell, probe concentration, and preferred scattering angles. In DWS 
and QELSS, the intensity autocorrelation function is obtained with a conventional digital 
autocorrelator; the time scales instrumentally accessible to QELSS and DWS are determined 
by the same instruments and are the same. The fluids and probes used by DWS are the 
same as those studied by Hallett, et al.-^ and many successors. OPD and DWS differ 
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secondarily in the opacity of their samples and primarily in the assumptions and formalism 
used to interpret their spectra. The cognomen Diffusing Wave Spectroscopy refers to a 
particular model for interpreting the field correlation function of the multiple-scattered light. 
In standard treatments of DWS in diffusing systems, it is proposed that DWS determines 
the mean-square diffusive displacement X{t) 2 of individual particles during time interval t. 
DWS has also been applied to nondiffusive systems as described below. 

This paper treats the contribution of particle motion and scattering path fluctuations to 
DWS spectra. It is shown below that fluctuations X 2n , n > 1, in the particle displacement 
in principle make non-negligible contributions to DWS spectra. Furthermore, even for paths 
of fixed length, these fluctuations are mixed with fluctuations from path to path in the 
number of scattering events and in the total-square scattering vector. Fluctuations lead to 
interesting challenges for the interpretation of DWS spectra in terms of particle motion. 

II. PHYSICAL BASIS OF THE DIFFUSING WAVE SPECTRUM 

We first examine general consideration o 19 i 20 i 21 that permit calculation of the diffusing 
wave spectrum. In essence, in a DWS experiment the light is scattered along many different 
paths before it reaches the detector. The scattered field incident on the photodetector is 
the amplitude-weighted coherent sum of the electric fields of the light travelling over all 
multiply-scattered paths. To describe the light scattered along a single path P having N 
scattering events, we take the positions of the iV scattering particles to be rj(£) for % e (1, N). 
The source is located at r ; the detector is located at r^+i- Unlike the particle positions, the 
source and detector positions are independent of time. The light initially has wavevector ko; 
the light scattered from particle i to particle i + 1 has wavevector kj(i). The light emerges 
from the system and proceeds to the detector with wavevector k^- The wavevectors k 
and kN are the same for every scattering path and at all times. The k; for 1 < % < N are 
functions of time because they link pairs of particles, and the particles move. Paths differ 
in the number of scattering particles and the ordered list of particles by which the light was 
scattered. For the scattering along a specific path P the phase shift is 

N 

A P <P{t) = -k„(t) ■ r - J2 r ^) ' ( " k i-iW ) + r ^+iW ' M*) ( 3 ) 

i=l 
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and the total field scattered along the available paths is 



E(t) = ^£ P exp(iA P 0(t)). 



(4) 



The summation is taken over the list of all allowed paths. The list of all allowed paths is the 
set of all ordered lists, of any length greater than zero, of particles in the system, subject to 
the constraint that no two adjoining elements in a list may refer to the same particle. Ep is 
the scattering amplitude associated with a particular path P. The phase shift A P (p(t) and 
scattering amplitude depend on time because the particles move. Because of translational 
invariance, the distribution of exp(iAp0(t)) is flat modulo 2ir, and the distribution of Ep is 
well-behaved. E(t) is therefore the sum of a large number of independent nearly-identically 
distributed random variables, so it is reasonable to infer from the central limit theorem that 
E(t) has a gaussian distribution. Because the E(t) have joint Gaussian distributions, the 
intensity autocorrelation function S(t) is determined by the field correlation function of the 
multiply scattered light. 

This description of multiple scattering reduces to the standard description of single- 
scattered light if the number of particles along a path is constrained to N = 1. By comparison 
with that description, the r^t) are the optical centers of mass of the diffusing particles, 
the single-particle structure factors arising from internal interference within each particle 
being contained in Ep. The standard description of scattering in terms of particle positions 
continues to be correct if the particles are non-dilute, and is the basis for calculating QELSS 
spectra of concentrated particle suspensions. 

The field correlation function depends on the phase shift as 



Because the paths are independent, in the double sum on P and P' only terms with P = P' 
are on the average non- vanishing. For those terms, the change between and t in the phase 
is rewritten by applying: the definition eq HU of the phase shift; the requirement that the 
source and detector locations and the initial and final wavevectors are independent of time; 
the definition qj(t) = (t) — kj_i(t) of the scattering vector; and addition of in the form 
qi(0) • rj(t) — qj(0) • Tiit). Re-ordering terms gives 



g {1 \t) = (E(0)E*(t)) = ( J]£p(O)£p,(t)exp[iAp0(t) -iA P ,0(O)] 





W = ( E p(o)E P (t) ex P -i q»(o) • 6n{t) - s qi (t) ■ n(t) 
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where Sri(t) = Ti(t) — r.;(0) is the particle displacement, and 5qi(t) = cy,(t) — qi(0) is the 
change in scattering vector i between times and t. The distance between scattering events 
is generally very large compared to the distances over which particles move before g^(t) has 
relaxed to zero, and the | kj | always have the same magnitude, so q.; can only be altered by 
changing the angle between kj_i and kj. In consequence, it is asserted that 5qj(t) term is 
extremely small relative to the q«(0) term, because | <5q.;(t) | / | qj(0) | and | <5rj(t) | / | rj(t) | 
are of the same order. The scattering amplitudes within the Ep only change as the scattering 
angles in each scattering event change, so, by the same rationale, for each path Ep(t) is very 
nearly independent of time over the times of interest. Furthermore, while the total of the 
scattering vectors must match the initial and final wave vector, if the number of scattering 
events along a path is large the constraint on the total scattering vector has very little effect 
on the intermediate scattering vectors. 

Under these approximations, the field correlation function reduces to 

9 {1) (t) = (J2 I e p(0) T exp • &■<(*)) ) , (?) 



p \ i=i 

,19 



as shown by Weitz and Pine ■ . The outermost average (• • • ) extends over all particle posi- 
tions and subsequent displacements. No assumption has thusfar been made about interac- 
tions between scattering particles. 

In order to evaluate this form, Weitz and Pine^ impose three approximations. Essentially 
equivalent approximations are imposed by other reference o 20 ' 21 . Each approximation replaces 
a fluctuating quantity with an average value. The three approximations are: 

Approximation 1) The exponential over the sum of particle displacements can be factor- 
ized, namely 

( I E P (Q) I 2 exp ( -i <J<(0) ■ M*)J ) = E I ^(°) I' II < ex P (-^(°) ■ 5r *W ' 

P \ \ 1=1 /IP i 

(8) 

so that averages over particle displacements can be taken separately over each particle. 
The factorization is valid if the scattering points along each path are dilute, so that the 
particle displacements 5ri(t), Sr 2 (t), . . . are independent, because under this condition the 
distribution function P( N \5ri(t), 5r2(t), . . .) for the simultaneous displacements of the iV 
particles of a path factors into a product of iV single-particle displacement distributions 
P ( - 1 \5r i (t)), one for each particle. In representative DWS experiments, mean free paths 



for optical scattering are reported as hundreds of microns^!, while the effective range of 
the very-long-range interparticle hydrodynamic interactions is some modest multiple of the 
particle radius or a modest number of microns, so the typical distance between scattering 
events is indeed far larger than the range over which particles can influence each other's 
motions. 

The assertion that the scattering points along each path are dilute does not imply that 
the scatterers are dilute. If the scattering cross-section of the scatterers is not too large, a 
given photon will be scattered by a particle, pass through many intervening particles, and 
finally be scattered by another, distant particle. The scattering particles may themselves 
be concentrated. The physical requirement on dilution used in ref. Q is that the mean free 
path between serial scattering events of a given photon is much longer than the range of the 
interparticle interactions. The motion of each scatterer in a path may very well be influenced 
by its near neighbors, but those near neighbors are almost never parts of the same path. 
While two very nearby particles could be involved in one path, this possibility involves a 
small fraction of the entirety of paths. Indeed, if most paths P included pairs of particles 
that were close enough to each other to interact with each other, the above description of 
DWS spectra would fail completely: Eq[7| would not follow from eq El because the terms of 
eqlHlin Sqi(t) would become significant. 



Finally, refs. 
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propose that the average over the individual particle displacements 



may be approximated as 

(exp (-icu(O) • 6n(t))) « exp {-c^XW/2) . (9) 



where X(t) 2 = ((g, • 5ri(t)) 2 ). A variety of rationales for eq El appear in the literature, as 
discussed below. Pine, et aL 21 show how eqElcan be replaced for nondiffusing particles when 
an alternative form is known a priori. 

Approximation 2) Each scattering event has its own scattering angle and corresponding 
scattering vector | q» |, which are approximated via (exp(— q^X 2 /2)) — ► exp(— q 2 X 2 /2), so 
that the scattering vector of each scattering event is replaced by a weighted average q 2 of all 
scattering vectors. Note that each pair q;, q^+i of scattering vectors shares a common kj, so 
that (qj-q i+ i) ^ 0. So long as the particle displacements 5rj(i) and 5r i+ i(t) are independent 
the cross-correlations in the scattering vectors do not affect the calculation. 

Approximation 3) The number iV of scattering events in a phase factor 
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exp ^— i ^2f =1 qi(0) • 5ri(t) j is approximated as being entirely determined by the opacity 
of the medium and the length of each path, so that all paths of a given length have exactly 
the same number of scattering events. Weitz and Pino^ propose that light propagation is 
effectively diffusive, there exists a mean distance /* over which the direction of light propa- 
gation decorrelates, the number of scattering events for a path of length s is always exactly 
N = s/£*, and the distribution of path lengths P{s) between the entrance and exit windows 
of the scattering cell can be obtained from a diffusive first crossing problem. 

These approximations were^ combined to predict that the field correlation function for 
DWS is 

POO 

0dws(*) « / dsP(s)ex V (-klxWs/t), (10) 
Jo 

where /cnis the wavevector of the original incident light and a mean scattering angle is linked 



by Ref. 



19| to £*. The model was solved for a suspension of identical particles performing 



simple Brownian motion, for which 



X(t) 2 = 2Dt. 



(11) 



For a light ray entering a cell at x — 0, travelling diffusively through the cell to x — L, 
and emerging for the first time into the region x > L, ref. Q finds the distribution of path 
lengths, based on the diffusion equation with appropriate boundary conditions. There are 
different solutions depending on whether one is uniformly illuminating the laser entrance 
window, is supplying a point source of light, or is supplying a narrow beam of light that 
has a Gaussian intensity profile. For example, for uniform entrance window illumination eq 



16.39a of ref 

fi'bws W 



19l gives for the DWS spectrum 



L/t + 4/3 
z /£* + 2/3 



sinh 



Zo_f&t_ 



0.5 



0.5 



- ( — ) cosh 
3 \r 1 



0.5 




L /6t 

I* \~T 



0.5 



+ 



4 fQt 



3 V r 



0.5 



cosh 



L f6t 

T 



0.5 



(12) 



where z Q w I* is the distance into the cell at which light motion has become diffusive, and 



r = (Dkl 



2\-l 



(13) 



is a mean diffusion time. 
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Even for an underlying simple exponential relaxation, the fluctuation in the total path 
length from path to path has caused the field correlation function for DWS to be quite 
complicated. If eq M were correct, then-ignoring the other approximations noted above- 
inversion of eqEDat a given t could formally obtain from flfowsW a va l ue f° r r a ^ that t. 
From r, eq IT31 formally gives a time- dependent D, whose frequency transform converts that 
time-dependent D into the frequency- dependent storage and loss moduli. 

III. FLUCTUATION CORRECTIONS 

Each of the above three approximations replaced the average of a function of a quantity 
with the function of the average of that quantity. Such replacements neglect the fluctuation 
in the quantity around its average, which is harmless if and only if the function is purely 
linear in the quantity being averaged. However, g^wstf) is an exponential of the fluctu- 
ating quantities. It is inobvious that fluctuations in its arguments can be neglected. We 
first consider the separate fluctuations in X, q 2 , and N and then demonstrate their joint 
contribution to a DWS spectrum. 

A. Fluctuations in Particle Displacement 

To demonstrate the relationship between the single-scattering field correlation function 



On the rhs, the average of the sum is the sum of the averages of the individual terms. By 
reflection symmetry, averages over terms odd in Ti(t) vanish. In the even terms, components 
of r-j that are orthogonal to q are killed by the scalar product. Without loss of generality, 
the x-coordinate may locally be set so that the surviving component of the displacement 
lies along the q axis, so q • r $ = qxi, leading to 



g^(q,t) = (exp(-iq-r i (t))) 



(14) 



and the mean particle displacements X(t) 2n , consider the Taylor series 




(15) 




oo 



(16) 



m=0 
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with (x™) = X n . The series on the rhs of eq HB] is the series ^^1 q 2n M n /n\ in q 2 , with 
expansion coefficients 

M = 1, (17) 



M x = — , 
2 ' 



(18) 



M 2 



X 4 
12' 



(19) 



X 6 
120' 



(20) 



etc. 



Equation Uni may also be written as a cumulant expansion exp(J^ n (— q 2 ) n K n /n\), the 
being cumulants. Expansion of the cumulant series as a power series in q 2 and comparison 
term by term shows 



9? fa*) = ex P 



2 X 2 4 (X 4 -3X 2 

*v - g - 



24 



5 



; (30X 2 - 15X 2 X 4 + X 6 ) 
720 



(21) 



with K 1 = ^ 1 K [ 



(X 4 -3X 2 ) 



2 • 
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, etc., the K n and X 2n being time-dependent. The cumulants 
K n here differ from the cumulants in Section D, below. If the distribution function for X(t) 
were a Gaussian, then all cumulants above the first would vanish (e.g., K 2 = 0), and eq l2*T1 
would reduce to eqs|21and|ni However, X(t) almost never has a Gaussian distribution in 
real systems. 



B. Fluctuations in the Scattering Vector 

For a single scattering event in which the light is deflected through an angle 8, the 
magnitude of the scattering vector is 

q = 2k sm(9/2), (22) 

where k = 2irn/\, with n the index of refraction and A the light wavelength in vacuo. 
Cumulant expansions of spectra depend on powers of q 2 as 

q 2 = 2k 2 (l -cos(0)). (23) 
11 



For DWS, the q 2 at the scattering points are independent from each other. All scattering an- 
gles are permitted. The average over all scattering angles comes from an intensity-weighted 
average over all scattering directions. For larger particles, the scattering is weighted by a 
particle form factor but always encompasses a non-zero range of angles. For small particles, 
no 9 is preferred. However, scattering from small particles is not isotropic, because light is 
a vector field. Scattering from small particles is described by dipole radiation, whose am- 
plitude is proportional to sin(^), ip being the angle between the direction of the scattered 
light and the direction of the polarization of the incident light. 

While it is true that light emerging from a turbid medium is depolarized, turbidity 
depolarization reflects the presence of many different paths, each of which rotates incident 
linearly polarized light through a different angle. If linearly polarized light is scattered by a 
typical small particle or a dielectric sphere (the typical probe) of any size, the scattered light 
from that one scattering event remains linearly polarized, though with a new polarization 
vector. In a typical QELSS experiment the incident light is vertically polarized, the first 
scattering event ip is measured from the perpendicular to the scattering plane, and the 
polarization remains perpendicular to the scattering plane, so sin(^) = 1. Because in 
multiple scattering the scattering paths are not confined to a plane perpendicular to the 
incident polarization axis, in general sin('0) 7^ 1. The factor sm(ip) arises already in QELSS 
experiments in which the incident laser polarization lies in the scattering plane. For example, 
for HH scattering (as opposed to the common VV experiment) from optically isotropic 
spheres, at 9 = 90° the scattering intensity is zero. 

For small particles, a simple geometric construction relates 9 and ip. Namely, without 
loss of generality we may take the scattering event to be at the origin, the incident light to 
define the +x-direction with its polarization defining the z-axis, and the scattering vector to 
be in an arbitrary direction not confined to the xy-plane. Defining (9 S , <p s ) to be the polar 
angles of the scattering direction relative to z, the +£-axis lies at (9 X , 4> x ) = (vr/2, 0) and the 
angle addition rule gives cos(6 l ) = cos(6' s ) cos(6' z ) — sin(6' s ) sin(9 x ) cos(0< j — <p x ), so one has 

q 2 = 2£; o 2 (l-sin(# s )cos(0 s )). (24) 

The angle 9 S is not the scattering angle 9 of eq|2Sl For the mean-square scattering vector 
from a single scattering event in a DWS experiment in which all scattering angles are allowed, 
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g 2 = (An)- 1 I dn s sm(6 s )2k 2 (l - sm(9 s ) cos(0 s )) = -k%, (25) 

The next two moments are g 4 = ^-^o and q e = ^-k®. The corresponding cumulants in a q 2 
expansion are K\ = nk 2 l /2, K2 = — 3.083A;q and K 3 = 0.7473/cq. The second cumulant is not 
negligible with respect to the first, in the sense that the variance (| K2/K 2 |) 1//2 is ~ 1.12. For 
in-plane scattering in a QELSS experiment, the average of q 2 over all allowed scattering an- 
gles would not be given by eqESl It would instead be be proportional to J d9 sin(#) sin 2 (6'/2). 
The averages for QELSS and DWS differ because for QELSS only in-plane scattering arises, 
while in DWS out-of-plane scattering events are allowed and important, and because in 
QELSS with VV scattering the corresponding polarization weighting factor is unity, while 
in DWS the allowed scattering angles are polarization weighted by sin(0 s ). A calculation 
made in the inadequate scalar-wave approximation would overlook this distinction. 

C. Fluctuations in the Number of Scattering Events 

In the standard treatment of photon diffusion in a DWS scattering cell, the path length 
distribution is computed by envisioning photons as random walkers, and solving the diffu- 
sion equation as a first-crossing problem to determine the distribution of path lengths. A 
pathlength s is approximated as containing precisely s/£* steps, the fluctuation in the num- 
ber of steps arising entirely from differences in the lengths of the various paths. MacKintosh 
and John 2 ^ present an extended treatment for the path length distribution, using a diffusion 
picture and saddle point methods to establish the mean N number of scattering events and 
the mean-square fluctuation in that number, as averaged over all path lengths. They treat 
separately paths involving few scattering events, for which a simple diffusion picture does 
not accurately yield the distribution of scattering events. 

In addition to the fluctuations in N arising from fluctuations in s, for a path of given s 

there are also fluctuations in N that arise because £* is only the average length of a path. 

While a path of length s on the average contains N(s) = sjl* scattering events, for paths 

2 

of given physical length s there will also be a fluctuation ([8N(s)) ) = N 2 (s) — N(s) in the 
number N(s) of scattering events. Scattering is a rate process linear in path length, so it 
is governed by Poisson statistics. For paths of fixed length, the mean-square fluctuation is 
therefore linear in the number of events. 
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D. Joint Fluctuation Effect 



For fluctuations within an exponential of a multilinear form, the case here, the error due 
to neglecting the fluctuation is readily obtained. Namely, for a function 

f(a) = (exp(—ax)) x , (26) 

where (• • • ) x denotes an average over the distribution of x, a Taylor series expansion gives 

/(a) = f;tf* (27) 
i=o % - 

Here x l = (x l ) x is the i th moment of x. The function f(a) may equally be written 

f(a) = exp (fj *f\ . (28) 

The Ki are the cumulants, with Kq = 1 and K\ = x. The higher-order cumulants Ki, which 
are the coefficients of a 1 in a Taylor series expansion for exp[a(x — x)] in powers of a, give 
the effects of the fluctuation. If X ^ 0, 

K 2 =^2_x 2 , (29) 

and 

fr 3 = F - 3a^x + 2x 3 . (30) 

An alternative case in which X = 0, and /(a) is a series in a 2n , is solved as eqs lT4Tl2Tl 

Cumulant expansions are well behaved, and converge under much the same conditions 
that Taylor series expansions are convergent. A cumulant series is particularly interesting if 
f(a) is very nearly exponential in a, because under that condition the relaxation is driven 
by K\ and the higher-order Ki are often all small. Cumulant expansions have already been 
used implicitly in the above. For example, the form (exp(— zq ■ m 1 — 9 x 2 ^ is the 

lowest order approximant. 

Where do fluctuations (equivalently, higher-order cumulants) modify the field correlation 
function for diffusing wave spectroscopy? The variables with interesting fluctuations are 
the displacement X, the mean-square scattering vector q 2 , and the number N of scattering 
events in a scattering path. The quantity being averaged is (exp(— Nq 2 (AX(t) 2 )). Repeated 
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series expansions in N, q 2 , and (AX) 2 , through the second cumulant in each variable, using 



the methods of ref 



AdwsW = ( exp I -N 



16, lead to 



q 2 X(t) 2 (f{X{ty - 3X(t) 2 
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N 2 (q 4 



q 2 (N 2 



N 2 ) 



+ ... 



(31) 



Here N and N 2 are the average and mean-square number of scattering events for all paths. 
One could also take N and N 2 to refer to paths of fixed length s, with (■ • • ) s including an 
average over the path length distribution. The above equation may be contrasted with the 
form <7dws(*) = (exp(— iV(s) q 2 X(t) 2 )) s , eq. EH obtained by approximating N(s) = s/£*, 
N(s) 2 - N(sf = 0, and X(t) 4 - 3X(t)2 2 = 0. 

Equation l3~Tl shows only the opening terms of series in the fluctuations in X 2 , q 2 , and N. 
The first line of eq 30 reflects particle displacements as captured by an individual single- 



scattering event and iterated N(s) times. The second line reflects the fluctuations from 
path to path in the total-square scattering vector and the number of scattering events, the 



fluctuations being g 4 — g 2 " and N(s) 2 — N(s) . The time dependence of gj}ws(^) in the above 



arises from the time dependences of X(t) 2 and X(t) 4 . The term X 4 — 3X 2 (and terms not 
displayed of higher order in X) reflect the deviation of distribution of particle displacements 
from a Gaussian. If q 2 and N were non-fluctuating, the second line of eq 30 would vanish. 
Because q 2 and N do fluctuate, S^ws^) § arns additional time-dependent terms, not seen in 
eqslTUlandl2n but appearing as the second line of eq l3*Tl 

On the rhs of the lead term of the exponential is the shorter-time approximant 



exp(— (N) (q 2 )X(t) 2 ) of Maret and Wol^. The basis of the approximation is that the dis- 
tribution of N has a peak location (approximated as N), paths having approximately N 
scattering events dominating the distribution, thereby giving the lead term to eq|^ The 
lead term does not give the true initial slope of g^ s (t). Because the decay rate of a path 
increases with increasing N, the minority of paths having particularly large N are respon- 
sible for the initial slope of gjyvysW- The minority of paths having particularly small N are 
substantially responsible for the slow decay of g^ s (t) at long times. Nonetheless, as shown 
in the following section, there is a sense in which DWS is sensitive to small-displacement 
particle motions at early times. 
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IV. EFFECT OF PARTICLE POLYDISPERSITY 



As a concrete example of fluctuation effects, we treat the DWS spectrum of a bidisperse 
system containing two sizes of probe. The probes perform Langevin-model diffusion in a 
simple Newtonian solvent. In this very special case, the single-particle field correlation 
function becomes 

g {1 \q,t) = A 1 exp(-D 1 q 2 t) + A 2 ex V (-D 2 q 2 t). (32) 

Here A\ and A 2 are the scattering cross-sections for the species 1 and 2 particles, respectively. 
D\ and D 2 are the respective diffusion coefficients. We transform eq EH into the canonical 
form of eql^by taking the exponential of the Taylor series of the logarithm of eqE21 namely 



g<H q ,t) = (A, + A 2 ) exp { - — tq + ^ + ^ t q + . . . j . (33) 

The mean-square displacement is determined by the average diffusion coefficient 

— {A X D X + A 2 D 2 ) . , 

(Ax + A 2 ) 

The second cumulant of the squared-displacement-distribution is determined by mean-square 
range of diffusion coefficients 

(AD 2 ) = AM A Dl ~ D f. (35) 
V 1 2{A l + A 2 ) 2 V ; 



By comparison with eqEH one finds X(t) 2 /2 = Dt and wtJ ~^ AW J = (AD 2 )t 2 . 

In interpreting g^(q,t), the average diffusion coefficient gives the initial linear slope of 
the spectrum, while on a semilog plot of g^ 1 ' against t the (AD 2 ) term gives a curvature of 
the spectrum away from its linear slope. It has been known since the earliest days of the 
QELSS technique that g^\q,t) is quite insensitive to weak polydispersity. Particles whose 
diffusion coefficients differ by a factor of two or three only lead to separable modes if the 
signal-to-noise ratio of the spectrum is extremely high, say, greater than 1000, permitting 
the spectrum to be followed through three orders of magnitude or more of relaxation. 

For a bidisperse probe mixture the field correlation function for a DWS spectrum becomes 

Sows W = (exp (-JV \f Dt + (AD 2 )tY/2 + ... 
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+ k 2 (g 4 - q 22 )D 2 t 2 /2 + (D q 2 t) 2 (N 2 - N 2 )/2] + ... 



) 



> 



(36) 



As an example of the implications of eqEEl consider a model bidisperse system similar to 
the system studied with DWS by Pine, et al.— , which contained mixtures of 198 and 605 nm 
polystyrene spheres in various concentration ratios. In convenient reduced units, the model 
spheres are given diffusion rates D\q 2 = 1 and D 2 q 2 = 3. To simplify the model, the two 
species contribute equally to scattering so A\ = A 2 , leading to D q 2 = 2 and (AD 2 )q 4 = 0.5. 
Pine, et al.— do not supply N for their experiment, but for cells of dimensions millimeters 
and mean free paths £* of hundreds of microns an interesting representative number might 
be N = 20. If the decay of the field correlation function can be observed over two orders 
of magnitude, then approximately N D q 2 t ~ 5, implying t « 0.13 in reduced units at 
the largest t observed. For this t, one has (AD 2 )qH 2 /2 < 0.01. Over the short range 
of times covered by the experiment, the deviations from single-exponential behavior at a 
signal-to-noise ratio of 100 is unobservably small . 

In this particular model, the time at which the spectrum has decayed to virtually to zero 
is so short that the second spectral cumulant (AD 2 ) has not yet contributed measurably 
to the spectral relaxation. The width of the path distribution leads to deviations from a 
simple exponential at very long times (paths with small N) and short times (paths with 
large N), but over the narrow range of intermediate times at which the DWS spectrum 
can be obtained, the spectrum is dominated by the first cumulant D q 2 . In agreement 
with these considerations, ref0 reports their mixture spectra are consistent with a single 
relaxation time. This result corresponds to the well-known QELSS result that mixtures of 
spheres of similar size have QELSS spectra that are very nearly pure exponentials. If the 
distribution of diffusion coefficients were made adequately wide, the terms in (AD 2 ) would 
become observable in the DWS spectrum. However, the terms in (AD 2 ) are determined 
by the mean-fourth power of the particle displacement, so the appearance of these terms 
in the observable spectrum would mean that the spectrum has ceased to determine the 
mean-square particle displacement. 

If the solution were polydisperse rather than bidisperse, the X(t) 2n would be more elab- 
orate. However, an essentially arbitrary relaxation spectrum can be written 




(37) 
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where A(T) is the normalized amplitude for relaxation at T, and as shown by KoppelM g(t) 
has a cumulant expansion 



where K\ = (T) = J dTA(T)T is the time-independent intensity-weighted average relaxation 
rate. Under circumstances numerically similar to the model, the DWS spectrum would 
only be sensitive to K±, because the DWS spectrum decays essentially to zero before K 2 = 
(r 2 ) — (r) 2 perturbs g{t). This rapid decay is the dark side of the optical level arm N 
advantage; just as the DWS spectrum is sensitive to very small particle motions, so also it 
does not readily see motions over larger distances. In contrast, the QELSS spectrum persists 
out to considerably longer (in natural units (D g 2 ) -1 ) times, so the second cumulant is more 
readily seen in QELSS spectra. If K 2 were substantially larger, it could perturb the DWS 
spectrum before the DWS spectrum relaxed. 

V. EFFECT OF PARTICLE INTERACTIONS AND OTHER PARTICLE MO- 
TION EFFECTS 

As seen above, under reasonable approximations, the phase shift for multiple scattering 
along a single path factors into a product of single-particle terms 



each of which refers to the motion of a single particle. The factorization is permitted because 
the mean path between scattering events is much larger than the distance over which particle 
motions are correlated, so that (Sri(t) ■ <5r i±1 (t)) = 0. While the scattering particles are far 
apart from each other, their physical neighbors that are not in the same scattering path may 
perturb their diffusive motions. In understanding the effects of interparticle interactions on 
QELSS and DWS spectra, it is important to recall 2 ^ that two physically distinct diffusion 
coefficients microscopically characterize the motion of diffusing macromolecules in non-dilute 
solutions, as studied with light scattering. One of these, the mutual diffusion coefficient, 
characterizes the relative motion of pairs of diffusing particles and describes the diffusion of 
particles down a macroscopic concentration gradient. The other of these, the single-particle 
or self diffusion coefficient, characterizes the motion of single particles through a uniform 
background. 




(38) 




(39) 
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Comparison with the modern literature on diffusing Brownian particle a 24 ^ 25 makes clear 
that QELSS is routinely applied in two experimental modes that correspond to the two 
diffusion coefficients. First, for single scattering from a solution of particles, all of which 
contribute equally to the scattering, one has for the field correlation function 

/mm \ 

<4%*) = ( E E ex P Hi ' (*<(*) " r ^°))) ) ( 4 °) 
\ i=i j=i I 

Here i and j independently label the M particles in the system. In this mode, QELSS 
determines the dynamic structure factor gm(q,t) and the mutual diffusion coefficient D m of 
the diffusing particles. D m follows from the first cumulant, namely 

n 2 ,. dg£\q,t) 

D m q =-hm — . (41) 

t->o dt 

Second, for scattering from a dilute solution of scattering particles in a perhaps- 
concentrated solution of non-scattering particles, QELSS determines to good approximation 
the single-particle correlation function 

gW(q,t) = M- 1 ^E ex P Hq " (n(t) - Ti(0)))\ . (42) 

Equation 021 is also the physical basis for optical probe diffusion. The single-particle cor- 
relation function that determines the DWS spectrum differs from eq 021 only in that q is 
different for each scattering event. From g^\q, t) the self diffusion coefficient of the diffusing 
particles follows as 

n 2 r dgj l) (q,t) 

D ° q = -^—dt~- (43) 
In nondilute solutions, the diffusion coefficients D s and D m are both modified by the 
direct and hydrodynamic interparticle interactions, but not in the same way. Interparticle 
interactions also contribute to the higher time cumulants of gm (q,t) and gi 1 \q,t), so for 
nondilute Brownian particles QELSS spectra are not exponential in time 2 ^. Correspondingly, 
for nondilute Brownian particles the displacement distribution function is not a Gaussian, 
because particle displacements are non-random: Successive displacements of Brownian par- 
ticles are correlated with each other because interparticle forces have long correlation times. 
The assertion that Brownian particles in complex fluids have non-Gaussian displacement 
distributions has been unambiguously confirmed by the direct experimental measurements 
of Apgar, Tseng, et a/ . 26 ' 27 
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The mutual and self diffusion coefficients are usefully written as averages over the hydro- 
dynamic interaction tensor hu, which describes the retardation in the motion of particle i 
due to the presence of neighboring particle /, and the interaction tensor T y -, which describes 
the motion induced in particle i by a force on particle j. For the drift velocity of particle 
i due to forces on particles j, one has^ 

N 

v Di (t) = J2^ F r ( 44 ) 
i=i 

Here the sum is over all N diffusing particles, /iy is the mobility tensor for the ij pair, and 
is the force on particle j. The mobility tensors are related to the hydrodynamic interaction 
tensors by 



1(1 + $>, + ...) (45) 

and 



/'-, -HT,, • ...) (46) 

Jo 

Here f a is the single particle drag coefficient, I is the identity tensor, and in eqEBH ^ j- 
For spherical particles, the hydrodynamic interaction tensors can be written as power series 
in a/rij, where a is a sphere radius and r y - is the distance between particles % and j. In 
particular, 

T« = ~[I + fyf«] (47) 

is the Oseen tensor approximation to T, with r y - being the unit vector pointing from particle 
i to particle j, while the self-term approximation corresponding to the Oseen tensor is 

, )' 

4 w, 



hu = -V ( — ) ( 48 ) 



Averages over the hydrodynamic interaction tensors, with due attention to the Brownian 
motion of the particles, fluctuation-dissipation requirements, and the correct interpretation 
of the nominal short-time limit of the time derivatives defining the diffusion coefficients 



givers 



+ exp(iq ■ r i£ )iqV e : (b u + T ie ))) (49) 
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and 

D s q 2 = (^f(q 2 + J2 (q-W-q))\. (50) 
\ Io i+i=\ I 

In the above equations, S(k, 0) is the static structure factor, is taken with respect to the 
coordinates of particle £, and the average is over all possible initial conditions. 

To which of these diffusion coefficients is diffusing wave spectroscopy sensitive? DWS 
responds to the correlation function of eq EH which is unmistakeably the single-particle 
correlation function seen in eq021 Correspondingly, DWS is sensitive to the self diffusion 
coefficient D s . The identification here that DWS measures D s has already been confirmed 
experimentally: Fraden and Maret^ and Qiu, et al.— used DWS to measure the diffusion 
coefficient D of polystyrene latex spheres as a function of sphere concentration, finding 
that the concentration dependence of the D measured with DWS agrees with theoretical 
expectations^ for D s . 

The above treatment of interaction effects differs from some prior analysis. MacKintosh 
and John^, their Section III, claim that in nondilute solutions one should replace eq I37H and 
their eq 3.6 

oo / n „ \ 

(E(t)E*(0)) oc J2 (II< / e'(r J ,t)e / (0,0)) ens exp(^ • Tj ) ) (51) 

n=l \ 1 Jt * I q 

with the rhs of their eq 3.7, 

N I b ^ I' (^5]«,/3exp(iq- [x Q (t) -x^(O)])) . (52) 



In eq|^ the j label n of the scattering particles, and time and space translational invariance 
have been used to start each particle at the origin at time 0, so rj is the displacement of j 
over time interval t. EqE21 describes a scattering event along some multiple scattering path. 
b(q) is a scattering cross-section, the sum is over all particles in the system, and x Q (t) and 
x j g(0) are the locations of particles a and (3 at times t and 0. The sum of particle positions 
in eqE21is the dynamic structure factor S(q,t). If this replacement were correct, which it is 
not, then DWS spectra would measure D m . 

However, in eq|^ comparison is only made between the position of the same particle at 
two times. In contrast, eq l5l?l terms involving the space-time displacements of distinct pairs 
of particles appear. A form like eqE21does appear in QELSS theory, in which scattering is 
coherent, so that the phase relationship between light rays scattered by a and (3 through q is 
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determined by the particle positions. Contrary to eqE21 in DWS, the paths leading from the 
laser to a and (3 are of independent, fluctuating length, so the fields scattered from particles 
a and (3 through q have independent random phases and can not interfere. Furthermore, 
in DWS a given particle scatters through q only if the previous and next particle along the 
scattering path lie along the pair of rays radiating from % that generate q. This constraint 
is far stronger than the constraint in QELSS, in which each particle scatters light in every 
direction consistent with the direction of the incident light. In a liquid, the DWS condition 
is generally not satisfied, so that even if one considers every multiple scattering path, only 
some particles scatter through any particular q. Contrarily, in eqE21 every particle in the 
system is assumed to scatter light to some other particle through each scattering vector q. 
Therefore, eq E21 is not a correct replacement for eq EHl DWS of non-dilute particles is in 
the first cumulant sensitive to the D s and not D m , in agreement with experimental results 
of Fraden, et al.— and Qiu, et al.— . 

Note that eq EH] does not require that particle motion be diffusive, only that the particle 
dynamics be known. Pine, et al.— discuss experiments on diffusing particles in shear flow, in 
which the mean-square particle displacement leading to the intermediate-time relaxation of 
the DWS spectrum includes both diffusive motions (which scale as q 2 t l ) and ballistic shear 
motions (which scale as q 2 t 2 ). Because the q and t dependences of these motions are both 
known a priori, the mathematical processes used to reach eqEHlfor diffusive particles lead 
equally to the DWS spectrum found experimentally for diffusing sheared particles, as shown 
by Pine, et al.— . 

MacKintosh and John^ proposes that eqE2]can be written 

gi 1 \q,t)^ex P (-q 2 W(t)) (53) 

with W(t) being the time-dependent mean-square particle displacement. They^fi cite Hess 
and Klein^i for this approximation, which is mathematically inconsistent with the exact eqliHI 
in that it is missing the terms in q 2n , n > 1. Hess and Klein^i discuss in detail the differences 
between exp(— q 2 W(t)) and gi^iq^t), carefully emphasizing that they treat W(t) instead of 
treating "... the full self-diffusion propagator, which is a complicated function of space and 
time. . . " . Hess and Klein further observe that gi {q, t) is the generating function not only 
for W(t) but also for the higher-order moments of the particle displacement distribution, 
those being the higher moments seen in eq[^ The use of W(t) may be traced back further 
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to Boon and Yip^ 2 - who give an expansion equivalent to eq |^ and note that even though 
g^\q, t) is not rigorously a Gaussian "it may be a good approximation to treat it as the 
sum of a Gaussian and a correction term" . The literature therefore does not support eq 
as a correct representation of eq WR 

In the discussion above on scattering from polydisperse systems, it was shown that DWS 
is typically sensitive only to the first cumulant of gi(q,t). At larger times, the higher 
cumulants of g^\q, t) are important, but before those times are reached the DWS spectrum 
may decay to zero. If the second and higher spectral time cumulants are sufficiently large, 
the deviation of the spectrum from a single exponential would visible. However, the higher 
time cumulants correspond to the higher moments of the displacement distribution function. 
If the single-particle function is not a pure exponential, then it does not reflect the mean- 
square particle displacement. 

VI. ANALYSIS 

In this paper, we treat the time dependence of diffusing-wave spectra. We demonstrate 
that the time dependence of DWS spectra arises, not only from the mean-square particle 

2 

displacement X(t) 2 , but also from deviations X(t) 4 — 3X(t) 2 from a Gaussian displace- 

2 

ment distribution, and also from higher powers X(t) 2 of the mean-square displacement. 
This result does not differ from the corresponding result for QELSS spectra, in which the 
time dependence of g^\t) arises not only from X(t) 2 but also from higher powers X(t) 2m 
of the mean-square displacement. Just as it is erroneous except as a crude approximant 
to write exp(— q 2 X{t) 2 ) for the general QELSS spectrum, so also it is erroneous to write 
exp(— (N) (q 2 )X(t) 2 ) for the general DWS spectrum. Furthermore, even though fluctuations 
in q 2 and N depend only slowly on time, the fluctuations couple to the strongly time- 
dependent X(t) 2n and thus to the time dependence of the spectrum at short times. The 
analysis of spectra of bidisperse systems shows that the approximant's error can be less 
serious in DWS spectra; namely, a DWS spectrum decays (N) times more rapidly that a 
QELSS spectrum, so a DWS spectrum may decay to zero before higher time cumulants 
become significant. 

The effect of fluctuations in N and q 2 on DWS spectra, phrased as deviations from eq 
IT2*| has been examined by Durian. Durian^ reported Monte Carlo simulations for photons 
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making random walks through a scattering slab. He computed the path length, number of 
scattering events, and sum of the squares of the scattering vectors for each path. These 
simulations determined fluctuations in the number of scattering events and the total-square 
scattering vectors, and determined the non-zero effect of these fluctuations on gjjws^)- 
Durian found that the fluctuation in the total square scattering vector Y increased more 
slowly than linearly with increasing pathlength s. A slower-than-linear increase is expected 
for a fluctuating quantity, and does not imply that the second cumulant of Y is negligible 
for long paths. From Durian's 33 simulations, an inversion of g^(t) via eq^Jto obtain X{t) 2 
has systematic errors, because eqH2]is inexact. Fluctuations described here and measured 
by Durian contribute significantly to the field correlation function. Durian demonstrates 
circumstances under which fluctuations in iV and q 2 only have effects of some small size on 
g^\t). He^ 3 . explains how Monte Carlo simulations could be used to overcome the effect 
of fluctuations, so as to make determinations of particle motion more accurate than those 
given by eq IT2l 

Discussions of light scattering spectra are sometimes referred to the Central Limit The- 
orem. The Central Limit Theorem provides that if a single random variable is constructed 
as a sum of a large number of identically distributed subsidiary random variables, then, as 
the number of subsidiary random variables becomes large, the distribution of the sum vari- 
able tends toward a monovariate Gaussian distribution. This theorem might be applied to 
describe the distribution of values of the field scattered by a large volume of solution at a 
single time, or the distribution of changes in the scattered field between any two times. 

However, the utility of the Central Limit Theorem is limited: 

First, the theorem requires that all subsidiary variables be identically distributed. In 
many cases of interest, different subsidiary variables in the sum have different distributions. 
Interesting cases in which the subsidiary variables are not identically distributed include: 
(1) The scattering Brownian particles have a bidisperse size distribution. In this case, 
the distribution of particle displacements is different for small and for large particles. (2) 
Different particles move in different environments. For example, the underlying complex 
fluid is approaching a critical point and has large long-lived local concentration fluctuations, 
so that particles moving in different fluctuations experience local media having different 
viscosities. (3) The underlying complex fluid is viscoelastic, so that the local viscoelastic 
properties of the fluid are in part determined by the local shear history, i.e., by how far each 
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particle has moved during prior times. For example, in a viscoelastic fluid, particles that 
had recently moved a greater-than-average distance might have perturbed the surrounding 
fluid more than their immobile neighbors would have, so at later times the resistance to their 
motions might differ from the average resistance, and the distribution of their displacements 
would differ from the average distribution. 

Second, the Central Limit Theorem gives the distribution of a single variable (which 
might be the difference between two other variables), but does not guarantee that the joint 
distribution of three random variables (for example, the values of the scattered field at three 
times), each pair of which has a joint Gaussian distribution, has an n-variate joint Gaussian 
distribution. Doob^ shows that to move from the Central Limit Theorem result for two 
random variables, to the results that three variables are jointly Gaussianly distributed, one 
needs an additional condition on the evolution of the subsidiary variables, namely that the 
subsidiary variables are described by a Markoff process. A simple example of a complex 
fluid system in which particle motions are not described by a Markoff process is example 
(3) of the previous paragraph. If one considers the displacements of particles between a trio 
of times t\ < £2 < t$, the displacement of particles between any two of those times might 
have a Gaussian distribution, but the three-fold distribution of particle displacements need 
not be a trivariate Gaussian, because the distribution of displacements between times £2 
and £3 might have a complicated dependence on the displacement between times £1 and £2. 
Note also the recent studies of Lemieux and Duria n 34 ^ 35 on intermittent dynamic processes 
leading to non-Gaussian scattering behaviors. 

Fortunately, there is an test that determines whether the Central Limit Theorem leads to- 
ward a calculation of the spectrum: It is a simple corollary^ of Doob's Theorem^ that if the 
particle displacements have identical Gaussian distributions, and successive displacements 
are all independent from each other, as is required for the particle motion to correspond to 
the Langevin model, then the QELSS spectrum is entirely determined by the mean-square 
particle displacement and is a pure exponential characterized by a time-independent diffu- 
sion coefficient and, correspondingly, a frequency-independent viscosity. Contrariwise, if the 
QELSS spectrum is not a pure exponential, the QELSS spectrum is not determined by the 
mean-square particle displacement. 

Even in the special case in which the QELSS spectrum is a simple exponential depending 
only on X(t) 2 , the DWS spectrum does not simplify. From eq 30, even if particle motions 
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are entirely characterized by X(t) 2 , so that X(t) 4 — 3X(t) 2 and similar higher-order terms 

2 

all vanish, the fluctuations in N and q 2 cause the DWS spectrum to depend on X(t) 2 
and higher order terms. In this special case that probe particles have identical Brownian 
displacement distributions, it might still be possible to extract D from a DWS spectrum. 
However, if the QELSS spectrum were not a simple exponential, for example because the 
optical probes were bidisperse, then the effect of the N and q 2 fluctuations is to make it far 
more challenging to extract the characteristics of the probe motion from a DWS spectrum 
than from a QELSS spectrum. The N and q 2 fluctuations mix fluctuations in the mean-ra th 

2 2 

displacements, such as X{t) 4 — X{t) 2 , with X 2 as seen in eq 30. 

VII. WHEN DOES DIFFUSING WAVE SPECTROSCOPY DETERMINE THE 
MEAN-SQUARE PARTICLE DISPLACEMENT? 

The underlying issue is interpreting the spectrum of light that has been scattered repeat- 
edly, perhaps many times, by a suspension of optical probe particles diffusing in a simple 
or complex fluid. The multiple-scattered spectrum is an elaborate average over sums of 
uncorrelated single-scattering events. The best that one can possibly do in interpreting a 
DWS spectrum is to make a perfect deconvolution of the averages over path length, number 
of scattering events, and scattering vectors. A perfect deconvolution would determine from 
the DWS field correlation function s^wsW ^ ne single-scattering field correlation function 
gs 1 \q,t). Even with a perfect deconvolution, the information from DWS spectra can be no 
better than the information in the single-scattering field correlation function gi 1 \q,t). We 
first consider what information is present in the single-scattering spectrum, and then con- 
sider additional issues that arise in attempting to deconvolve g^ s {t) to determine gs(q, t). 

A. Interpretation of gi 1 \q,t), However Obtained 

As seen from eqEU the single-scattering gi^iq^) is determined not only by X 2 but also 
by all higher moments X 2n , n > 1. The higher moments appear in a variety of combinations, 
the combinations being non-zero except in the special case that the displacement distribution 
P(X) is a Gaussian in X. In the special case of a Gaussian P(X), gi (q,t) reduces to a 
Gaussian in X 2 . Contrariwise, if P(X) is not a Gaussian, then the higher even moments of 
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A all contribute to gi 1 \q,t). 

P(X) can be measured directly via video microscopy, at least in systems in which par- 
ticles do not move too quickly. For example, Apgar, Tseng, and collaborator a 26 i 27 report on 
mesoscopic probe particles in water:glycerol and in aqueous actin, actimfascin, and actin: a- 
actinin mixtures. In this work, video light microscopy was used to make repeated mea- 
surements of the positions of large numbers of particles at many times; the distribution 
of particle displacements during various time intervals was computed. Probe particles in 
water:glycerol show a Gaussian displacement distribution. Probe particles in the protein 
solutions have markedly non-Gaussian displacement distributions, with a displacement dis- 
tribution far wider (in terms of X 2 /X 2 ) and more skewed than observed for particles in 
water:glycerol. The direct measurements of Apgar, et alj2& and Tseng, et al^L show that 
any assumption that P(X) for probes in water: actin is a Gaussian is quantitatively incor- 
rect. Correspondingly, X 2 (t) does not characterize probe motion in these systems. Analyses 
of DWS spectra of probe motion in water: actin systems based on eq. ^] are be invalid, 
because the underlying assumptions behind eq. ^Jare not satisfied. 

Equation El for g < a\q-,t) shows that P(X) and gji(q,t) are linked by a spatial Fourier 
transform between X and q. If P(X) is a Gaussian in A, then its transform gi i (g, t) must be 
the Gaussian exp(— q 2 X 2 ) in q that was assumed in the derivation of eq^J Contrariwise, 
if gi^iq^t) is not a Gaussian in g, then by the same token P(X) is not a Gaussian in 
A, and (^(g, depends not only on A 2 but on the higher moments of A. There is a 
considerable literature on optical probe diffusion as studied with QELSS, often by applying 
mode decomposition or related spectral analysis methods to gs\q, t). For probes in HPC 
solutions, Streletzky, et al .? 7 ^ show: While the relaxation rates of some spectral modes 
scale as q 2 , for other modes the relaxation rates are not linear in q 2 . From equation EJ 
the component of the displacement distribution P{X) corresponding to a mode that does 
not relax as exp(— aq 2 ) is necessarily not Gaussian in A. It would be incorrect to infer 
the viscoelastic properties of systems with relaxations whose relaxation rates do not scale 
linearly in g 2 by using eq El or El because particle motions in these systems would not satisfy 
the assumptions on which these equations are based. 

The relaxation of an arbitrary mode can formally be written exp(—tq 2 D(q,t)). With 
simple diffusion, D(q,t) is a constant having the trivial q and t dependences q°t°. If D(q,t) 
has nontrivial dependences on q and t, an average over q does not factorize as g 2 -D = q 2 D 
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because D is a function of q, contrary to the implicit assumption that the exponent on the 
rhs of eq M was purely quadratic in q. 

What are the special cases in which eq provides a correct description of the single- 
scattering field correlation function? Two are readily identified. First, the probes might be 
a monodisperse suspension that diffuses in accord with the Langevin equation, as described 
by Berne and Pecora^. Polystyrene latex spheres in water:glycerol afford an example. The 
solvent is a simple Newtonian fluid having no viscoelastic memory on observable time scales. 
In this case, the distribution of particle displacements is Gaussian, and 



with D a constant. The diagnostic for this case is that a plot of \og(g s (t)) is a straight 
line, starting at the smallest observable time and extending out until the signal fades into 
the noise. 

Second, the probes might be diffusing in a viscoelastic fluid that has identifiable longest 
time and distance scales on which relaxation occurs. Particle motions, over times and dis- 
tances much longer than the largest relaxation time and distance, satisfy the requirements 
of the Central Limit Theorem and Doob's Theorem. Over sufficiently large times and dis- 
tances, the probes perform simple Brownian motion. At long times \og(gi 1 \t)) against t 
becomes a straight line, from which a long-time X(t) 2 and D can be extracted. 

To extract a long-time limiting slope from \og(gi (i)), a long time linear limit must 
actually exist, a circumstance that is not guaranteed to arise. For example, careful QELSS 
studies show that in some polymer systems gs (t) decays at large t as a stretched exponential 
in time^. These systems have a continuous distribution of relaxation times but not a single 
longest relaxation time. In these systems, fitting a straight line to log(</i (t)) at large t is 
meaningless. 

B. Deconvolution of g-jjwa(q,t) to Determine gi (q,t) 

The reconstruction of the single-scattering field correlation function from the diffusing- 
wave field-correlation function faces a fundamental challenge. It is fundamentally impossible 
to reconstruct a general function of two variables q and t, namely gs (q, t) from a general 
function of one variable t, namely g^wsW' wnen t ne univariate function of t was generated 




(54) 
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via an average over q of the bivariate function. The issue is simple: information on the 
g-dependence of gj (q,t) is destroyed by the averaging process. This fundamental limit 
is mathematical, not physical, and arises from the information-theoretic consequences of 
taking the average over q. 

It might superficially appear that eq fT2l reconstructs gj (q,t) from fl^wsW- There is no 
reconstruction here. Rather, the g-dependence of gj (q,t) is taken to be known a priori to 
be ga(q,t) ~ exp(— q 2 Dt), for D independent of q. Equation IT2l then only needs to recon- 
structs the time dependence of gj 1 \q,t), and that only so far as the first time cumulant, 
which is possible. The q dependence of gj (q, t) does not need to be diffusive. Wu, et al.— 
treat reconstruction for diffusing particles in a shear flow, for which gs(q,t) has known 
diffusive and shear components. A g-dependence that is known a priori will support recon- 
struction. However, without a known form for the g-dependence of gs(q,t), reconstruction 
of gj (q, t) from gj} WS (£) be made. For probe particles in polymer solutions, a common topic 
of investigation, the careful work of Streletzky^ 7 -^ shows that the q dependence of gj 1 \q, t) 
is a complicated function of probe radius and polymer concentration that must be deter- 
mined by systematic experiments. Probe motion in nondilute polymer solutions is therefore 
fundamentally inaccessible to study by DWS as the method is presently constituted. 
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